Coherent control via interplay between driving field and two-body interaction in a double well 
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We investigate interplay between external field and interatomic interaction and its applications to coherent 
control of quantum tunneling for two repulsive bosons confined in a high-frequency driven double well. A 
full solution of the system is generated analytically as a coherent non-Floquet state by using the Floquet states 
as a set of complete bases. It is demonstrated that the photon resonance of interaction leads to translation of 
the Floquet level-crossing points, and the non-resonant interaction causes avoided crossing of partial levels. In 
the non-Floquet states, the bosons beyond the crossing points slowly vary their populations, and the resonant 
(non-resonant) interactions enhance (decrease) the tunneling rate of the paired particles. Three different kinds 
of the coherent destructions of tunneling (CDT) at the crossing, avoided-crossing and uncrossing points, and 
the corresponding stationary-like states, are illustrated. The analytical results are numerically confirmed and 
perfect agreements are found. Based on the results, an useful scheme of quantum tunneling switch between 
stationary-like states is presented. 
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I. INTRODUCTION 

Coherent control of quantum tunneling in a double well 
via periodical driving has been researched extensively from 
both theoretical and experimental sides |jl]-3]. Some interest- 
ing phenomena, such as the coherent destruction or construc- 
tion of tunneling (CDT or CCT) fllToll. chaos enhancing tun- 
neling 111 li I12II . and photon-assisted tunneling |13, 14] have 
been found. Many works focus on single- or many-particle 
systems. Few-particle systems are of a intermediate class be- 
tween the both systems, which deserve further study for com- 
prehensively understanding tunneling dynamics. However, in- 
vestigations on quantum control to two particles in a period- 
ically driven double well are extremely rare, except for the 
cases of two particles in a one-dimensional lattice which can 
be reduced to a two-site trap ifTsi IT6I1 or two particles in a 
driven double-well train IUTI - uOIi . Recently, some relevant re- 
searches have been completed for a non-diiven few-paiticle 
system ll2ll - l23ll . Several interesting phenomena of quantum 
tunneling were shown for two repulsive bosons in a non- 
driven double well, which include the Rabi oscillations and 
correlated pair tunneling ll24l - l26l . The interatomic interaction 
adjusted by the Feshbach- resonance technique [27] plays an 
important role in tunneling dynamics of the non-driven two- 
particle system. Here we are interested in the combined ef- 
fects of the driving and interaction on tunneling dynamics of 
double- well coupled two bosons. 

In this paper, we study coherent control of quantum tunnel- 
ing via the competition and cooperation between atomic inter- 
action and driving field for a pair of repulsive bosons confined 
in a periodically driven double well. In the high-frequency 
regime, we obtain a set of Floquet quasienergies and Floquet 
states of invariant population, which contain the quasi-NOON 
state, an interesting entanglement state ll28[l29ll . The general 
non-Floquet state of slowly varying population is generated as 
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a full solution which is a coherent superposition of the Floquet 
states. The quasienergies as function of the driving parameters 
are plotted for several different values of interaction strength. 
The quasienergy spectra exhibit that coniparing with the non- 
interacting case, the photon resonance fly\ leads to transla- 
tion of the level-crossing points, and the non-resonant interac- 
tions cause avoided crossing of partial levels. Subsequently, 
exploiting the coherent non-Floquet solutions, we investigate 
time evolution of the particle population and demonstrate that 
beyond the crossing points, the bosons slowly vary their popu- 
lations compared to the high-frequency driving. The resonant 
or non-resonant interactions can enhance or decrease the tun- 
neling rate of the paired particles. Three different kinds of 
CDT, respectively at the level-crossing points for the super- 
position of three states, at the avoided-crossing points for the 
superposition of two states and at the uncrossing points for 
the single Floquet states, are illustrated, which result in differ- 
ent stationary-like states of invariant populations. The analyt- 
ical results are confirmed by direct numerical simulations and 
good agreements are shown. As an application of the above 
results, an interesting scheme of quantum tunneling switch be- 
tween stationary-like states is presented by using CDT or CCT 
to close or open the quantum tunneling, which could be useful 
for the quantum control of two bosons in a double well. 

II. ANALYTICAL SOLUTIONS IN HIGH-FREQUENCY 
REGIME 

We consider two ultracold bosons confined in a periodically 
driven double well with the governing Hamiltonian 1 30, .3 1.1 

Hit) - ^ib\bi-blb2) + y(blb2 + blbi) 

+ ^[blbi(blbi-l)+blb2(blb2-l)], (1) 

where /?,(/?■) for i - 1,2 are the annihilation (creation) oper- 
ators in the /th localized state which may be a line superposi- 
tion of the lowest doublet of single-particle energy eigenstates 
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ll32ll . Their commutation relations read [bi,b^j] - 6ij. The pa- 
rameters y and U are the tunneUng coefficient between two 
wells and the interaction strength between two bosons. The 
function e(f) = sq cos(wf) describes the external ac field in 
which a> and eo are the driving frequency and amplitude, re- 
spectively. For simplicity, we have put H - I and taken the 
reference frequency wq = lO^i"' [33] so that eq, y, a> and U 
are in units of coq, and time t is normalized in units of Wq ' . 

Quantum state \i^(t)) of s^tem (1) can be expanded in Fock 
bases |0,2),|1, 1),|2,0) as Oj 

\m) = flo(f)|0,2) + fli(f)|l,l) + fl2(f)|2,0), (2) 

where \j} = \j,2 - - 0, 1,2) denote that j bosons re- 
side in left well, 2 - j bosons reside in right well; a jit) rep- 
resents probability amplitude of the system in y'-th Fock state 
which obey the normalization condition |flo(f)P + kiCOP + 
l«2(0P - 1- Inserting Eqs. (1) and (2) into the time-dependent 



Schrodinger Equation / 
pling equations 



— H{t)\ip{t)) results in the cou- 



;ao(0 = [f/ - e(f)]flo(0 + ^yai{t), 
iai(t) ^ ^yao(t)+ ^yaid), 
iaiit) ^[U + e(0]a2W + ^yai(t). 



(3) 



It is difficult to get exact analytical solutions of Eq. (3) 
with periodic driving s{t) = eocos(wf). However, we can 
construct the approximate analytical solution in the high- 
frequency limit, w » 1 and m, y <K w. To do this, we adopt the 
idea of reduced interaction strength ifisl to rewrite the interac- 
tion strength as [/ = no) + u for < m <K w, n = 0, 1 , 2, ... with 
u being the reduced interaction strength, and make the func- 
tion transformation ao(f) = boiOe' , ai{t) - bi(_t), 
aiit) = fe2(f)e"'^*''^"^'™*'", where bj(t)(j = 0, 1,2) are slowly- 
varying functions of time, then Eq. (3) becomes new equa- 
tions in terms of bj{t). Because of \bj{t)\ - \aj(t)\, the 
high-frequency limit implies that the occupied probabilities 
of Fock states slowly vary in time. Exploiting Fourier expan- 
sion exp[ J(+/e(f) - ma>)dt] - exp[+/(^) sin(wf) - incot] - 
Yj'^>__^J^n'(^)sxp[i{+n' - n)(L)t], we easily obtain the time- 
average of the rapidly oscillating function as jr±«(^) which 
are the +«-order bessel function. Under the high-frequency 
approximation, the rapidly oscillating function can be re- 
placed by its time-averaging value 1,30,1 such that the equations 
of b j(t) are transformed into 

ibo(t) = uboit) + J„{—)bi(t), 
a> 

ibiit) = j„{-)hit) + (-irj„(-)b2(t), 

CO CO 

Mt) = Ub2{t) + {-\yjn{-)bi{t). (4) 
CO 

In the calculations, we have employed the formula JT-ni^) - 
(-1)"^h(— ) for positive integer n, and the renormalized cou- 
pling coefficient 7„(^) - yf2yj„{^). Starting from Eq. (4), 
we obtain the interesting analytical solutions as follows. 



A. Quasienergies and Floquet states 

Because the time-dependent Hamiltonian (1) has the pe- 
riod T = we can make use of Floquet theory ll34ll to 
get its solution in the form |i/'(0) = e^'^'\(f(t)), where \(p{t)) 
is the Floquet state with the same period E is called the 
Floquet quasienergy which has been normalized in units of 
Hcoq. Noting the same period of the transformation func- 
tion e/(±'s(')-'"(<^)<i'' between a/t) and bj(t), to generate the Flo- 
quet states, we seek the stationary solutions [35] of Eq. (4) 
bo = Ae-'^', bi = Be-'^', bi = Ce"'^' with constants A, B, C 
obeying the normalization condition |A|^ + |Bp + |Cp = 1. In- 
serting these into Eq. (4), we establish the equations of the 
stationary solutions as 

{E - u)A - J„{—)B = 0, 

CO 

-U-)A + EB- (-1)"J„(-)C = 0, 

CO CO 



-(-l)"J„{—)B + (E- u)C = 0. 



(5) 



By solving Eq. (5), we obtain three Floqeut quasienergies £; 
and three sets of constants A/, Bi, Ci for / = 0, 1, 2. 

£o = «,Ao = ^,Bo = 0,Co = (-l)"+'^; (6) 
V2 V2 



1 2J„ 
El = -{u-k„),Ai^ 

2 V8-/„' + (« + kn)^ 



V872 + (m + fc„)2 



V8y2 + (M + ^„)2' 



(7) 



1 2J„ 

El = -(M-Hfc„),A2 = 

2 V8-/„' + (k„ - m)2 

B2 = -r.^r-n" 



+ {kn - m)2 



^l%Jl + (k„ - m)2 



(8) 



Here, the simplified parameter k„ - ^87,^(^) + u?- has been 

adopted, so that the Floquet quasiegergies are functions of 
driving parameters and interaction strength. Obviously, for 
M, y <K w we have Ei <K w. Thus as the functions of exp[-/£'/f] 
the stationary solutions bj(t) are slowly varying indeed. Given 
bj(t), the rapidly oscillating function aj(t) are obtained imme- 
diately. Substituting such cijif) into Eq. (2), Floquet states 
|i/';(f)) are expressed as 

|i/,,(f)) = e-'^''(A,e^^'"*"'>-'"'"|0,2)-HB,|l,l) 

H- C,e^^'"*"'^-'""'|2,0)) (9) 

for Z = 0, 1 , 2. Here n is a positive integer adjusted by the inter- 
action strength. If the system is prepared in the Floquet states 
initially, probabilities of the system in different Fock states are 
the constants |A;p, \Bj^ and |C/p, respectively, so that tunnel- 
ing of atoms between two wells is suppressed completely, i.e., 
CDT occurs. 
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B. Coherent non-Flouqet states 

Given the quasienergies and Floquet solutions of Eqs. (6- 
9), the principle of superposition of quantum mechanics in- 
dicates that the linear Schrodinger equation has the periodic 
or quasiperiodic non-Floquet solutions, which is coherent su- 
perposition of the Floquet states as a set of complete bases 
in a three-dimensional Hilbert space [Esl 1321). Noticing the 
relation between the rapidly oscillating probability amplitude 
aj{t) and the slowly varying function bj(t), the general non- 
Floquet state has the form 

= co|(/'o(0> + ci|iAi(f)> + c2|i/'2(0) 
= Uoim, 2) + a[(t)\l, 1) + fl^(0|2, 0) 

= b[^(t)e'^ sinM)-"™'|0, 2) + b[it)\l, 1) 



+ b'^{t)e- 



(10) 



where c/ (Z = 0, 1 , 2) are superposition coefficients determined 
by the initial conditions and normalization, and b'j(t) are the 
slowly varying functions 



b'^it) = 




-'^'" + ciAie 


+ C2A2e- 


-!&' 


(11) 


b\{t) = 


cqBqc 




-'^" + C2B2e- 




(12) 


b'^it) = 


coCoe 






-iE2t 


(13) 



with constants A;, Bi, Ci being given in Eqs. (6)-(8). Clearly, 
the functions a'j{t) and b'j{t) obey Eqs. (3) and (4), respec- 
tively, and have the same norm which is the corresponding 
occupied probability. Any b'j(t) is a superposition of three 
periodic functions with frequencies wj - Ei for I - 0, 1,2. 
When all the frequency ratios (jj'i/a>'i„(l + m) for any pair of 
I, m = 0,1,2 are rational numbers, any b'j{t) is a periodic func- 
tion, and any irrational frequency ratio implies that all the b'j{t) 
for 7 = 0, 1, 2 are the quasi-periodic functions ll35ll . 

The non-Floquet state of Eq. (10) with Eqs. (11-13) is a 
general solution with constants c/ being determined by the ini- 
tial conditions and normalization. As a example, we consider 
the two atoms reside in right well initially, which means the 
normalized initial state of the system to be \t//(0)} = |0, 2). In- 
serting this into Eq. (10) yields b[^(0) = l,b\(0) = b'^(0) = 0, 
then substituting these and Eqs. (6-8) into Eqs. (11-13) pro- 
duces 



m = -e- 



4kn 



4k„ 



b'^if) = (-l)«+iie-"'' + (-l)«^:i^e-^("-*^"^' 

2 



(15) 



(-ly 



+k„)t 



(16) 



Combining Eq. (10) with Eqs. (14-16), we arrive at the spe- 
cial non-Floquet state associated with the initial state \iff{0)} - 
|0, 2). Such a special state will be used, as an instance, to show 
the coherent control of quantum tunneling. 



III. EFFECTS OF DRIVING AND INTERACTION ON 
QUASIENERGY 

We have already obtained the Floquet quasienergies, Flo- 
quet states and non-Floquet states in high frequency regime. 
The general non-Floquet state of Eq. (10) is a coherent su- 
perposition of two or three Floquet states with two or three 
nonzero superposition constants c/ determined by the initial 
setup. It is worth noting that in Eqs. (11-13) the three 
quasienergies appear in the time-dependent phases of the 
probability amplitudes and directly affect the tunneling proba- 
bilities through the phase coherence. When the quasienergies 
of the Floquet states are different each other, the tunneling 
probabilities are periodic or quasiperiodic functions of time, 
which describe coherent population oscillation of the system. 
If all the three quasienergies are the same, Eqs. (11-13) mean 
invariant populations with constant \b'j{t)\, namely the level- 
crossing leads to CDT and stationary-like states. Similarly for 
the case £0 = £2 ^ Ei, the superposition state of two Flo- 
quet states |i/'o) and \t//2} describes the second kind of CDT, 
where the energy E-i avoids the level-crossing. Therefore, 
the analytical results render the the relation between level- 
crossing and CDT more transparent and the dependence of 
tunneling dynamics on the values of Floquet quasienergies 
more evident. From the analytical results of Eqs. (6-8) we 
know that the quasienergies are the functions of interaction 
strength and driving parameters. In this section, we study in- 
terplay between the external driving and interatomic interac- 
tion, through the quasienergy spectra. 

For a high frequency the expression U = nio+u with u <^ oj 
implies that « = denotes noninteracting or weakly interact- 
ing case, n > 1 correspond to strongly interacting case. Here, 
we will consider only the cases n = and n - \, due to the 
interaction U - to + u strong enough. From Eqs. (6-8) the 
Floquet quasienergies as functions of the driving parameters 
for several different values of interacting strength U and the 
same coupling coefficients y are shown as the circles in Fig. 
1 . We also calculate the quasienergy spectra numerically from 
Eq. (3) for the same parameters as those of the analytical re- 
sults. The numerical results are plotted as the solid curves of 
Fig. 1 and both the results consistently illustrate the following 
interesting phenomena. 

Photon resonance leads to translation of the level-crossing 
points. In Fig. 1(a) and Fig. 1 (c) we show that for the non- 
interacting and resonant interaction cases with u - Q and n - 
0, 1, quasienergy spectra exhibit exact level-crossings. Three 
quasienergies are same at the crossing points, Eq - E\ - E2, 
such that the probabilities \b'j(t)\^ (J = 0, 1,2) of Eqs. (11)- 
(13) do not change in time, and the CDT of the first kind oc- 
curs. In the absence of interaction (U = 0), the inset of Fig. 
1(a) shows the first level-crossing location xi = sq/lo = 2.405 
which is the first root of equation J'oix) - 0. While for the 
strong interaction (C/ = 50) the inset of Fig. 1(c) indicates 
the first level-crossing location yi = eo/w = 3.832 which 
is the first root of J'l (y) = 0. The results mean that photon 
resonance with U - cj results in translation of the location 
of level-crossing point from (sq/cj, U) - {x\, 0) to {yi, ai). 
This implies that in the strong interaction case the onset of 
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CDT requires a greater sq/w value, because of yi > xi. 

Non-resonant interaction causes avoided crossing of par- 
tial levels. In Fig. 1(b) and Fig. 1(d) we show that for u + Q 
and « = 0, 1 the crossing of three levels is replaced by that 
of two levels at the locations (so/ at, U) = (2.405, 2) and 
(3.832, 52), respectively, namely the avoided crossing of one 
level appears. Such avoided level-crossing points of partial 
levels correspond to driving parameters obeying J'niso/a)) - 
0, which are the points of closet approach of energies in quasi- 
energy spectra. At such points the coherence between phases 
Eot and Eit may exist in the non-Floquet state of Eq. (10) for 
the nonzero q, / = 0, 1, 2. However, for the constants ci = 
and Co, C2 9^ 0, the non-Floquet state becomes the superpo- 
sition of two Floquet states with the same quasienergy and 
time-dependent phase, so the probabilities \b'j{t)[^(j = 0, 1,2) 
of Eqs. (11-13) don't vary in time and the second kind CDT 
occurs. 
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FIG. 1 ! (Color online) Quasienergy spectra for the driving frequency to - 50, coupling 
coefficient y = 0.5 and interaction strength (a) U - 0; (b) C/ - « - 2; (c) (7 - - 50; 
(d) U = w + « = 52. Hereafter any quantity plotted in the figures is dimensionless, 
circles or circular points denote the analytical results and solid curves the numerical 
correspondences respectively, unless it is specially indicated. 



b'^iO) = l,b\(0) = b'^iO) = correspond to Eqs. (14-16) and 
other initial setups are associated with Eqs. (11-13). 



A. Manipulating tunneling rates 

At first, we consider the case in which the system parame- 
ters beyond level-crossing points and the bosonic population 
oscillates periodically. We study how to control the tunneling 
rates by setting and adjusting the values of interaction strength 
and driving parameters. 

Letting P jit) - \b'j{t)\^ be probability of the system in the 
j-th state, from Eqs. (14-16) we plot time evolutions of the 
probabilities for j - 0, 1,2, w = 50, sq - 2a) = 100 and 
y = 0.5, as shown in Fig. 2. In this figure we show that all 
the probabilities are slowly varying function of time compar- 
ing with the high-frequency driving. Defining tunneling time 
as the needed time of the paired bosons tunneling from state 
|0, 2) to |2,0), a rough estimate from curves 1, 2, 3 and 4 in 
Fig. 2(a) and the inset of Fig. 2(b) indicates non-resonant 
tunneling time of the paired particles for the weak interaction 
U = 2 to be about 125 and for the strong interaction U = 52 to 
be about 20, compared to 6 with resonant interaction [/ = 50 
and 14 without interaction. The results mean that photon res- 
onance iU - to) induces coherent construction of tunneling 
(CCT) which leads to increase of the tunneling rate of paired 
particles comparing with the non-resonance case. Especially, 
for the weak interaction (U - 2), Pi oscillates around nearly 
and its values are negligible. Thus the two bosons tunnel as 
pair between wells such that at any time, quantum state of the 
system is a superposition of states |0, 2) and |2, 0), which ar- 
rives at the NOON state with probabilities Pq - P2 - 1/2 
periodically. 

We also obtain numerical results from Eq. (3) with the same 
parameters and initial conditions fl()(0) = ^(')(0) = l,fli(0) = 
fl2(0) - ^j(O) - b'2{0) = as those of the analytical case, 
which are shown as solid lines of Fig. 2. The analytical results 
are in good agreement with the numerical results except the 
slight deviation for the case U = 52. 



IV. CONTROL AND SWITCH OF QUANTUM TUNNELING 

We have known that at the points without level-crossing 
the bosonic population oscillates periodically or quasi- 
periodically, and at the level-crossing points the initial pop- 
ulation can be kept. Now we investigate the coherent con- 
trol of quantum tunneling by applying the interplay between 
the external field and interatomic interaction. For the level- 
uncrossing case we study how to coherently manipulate the 
tunneUng rates by setting and adjusting the values of the sy- 
atem parameters. For the level-crossing and avoided cross- 
ing cases we perform control to the stationary-like states via 
CDT of different kinds. Then we apply these results to present 
a scheme for designing the quantum tunneling switch from 
a given state to different stationary-like states under CDT. 
Hereafter, all the analytical results are based on the gen- 
eral non-Floquet state of Eq. (10), and the initial conditions 



B. Preparing stationary-like states via CDT of different kinds 

When the tunneling rate is controlled to zero, CDT occurs 
and the stationary-like state of invariant population is pre- 
pared. Such stationary-like state may be a single Floquet state 
or the superposition state of three or two Floquet states, which 
correspond to CDT of three different kinds, respectively. 

CDT at the level-crossing points. At level-crossing point 
(si)/a>, U) = (2.405, 0) and for the high-frequency co = 50, 
time evolution of Po(0 = l^o(OI^ of Eq. (14) is shown as 
curve 1 in Fig. 3(a). It is observed from this curve that 
P(){t) maintains the initial value Po(Q) = 1 so that tunnel- 
ing between wells is suppressed completely and CDT of the 
first kind occurs. Under the same driving parameters, for the 
level-uncrossing points with strong interaction U - 50 and 
52, Po(t) still oscillates between and 1, as indicated by the 
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FIG. 2: (Color online) Time evolutions of (a) Poit), (b) Pi (/) and (c) ^2(0 for the parameters eq = 100, aj = 50 and y = 0.5, and U = (curve 1), t/ - 2 (curve 2), U - 50, (curve 
3) and f/ = 52 (cui've 4). The long-time evolutions of curve 2 in Figs. 2(a) and 2(c) are plotted in the inset of Fig. 2(b). 



curve 2 and curve 3 of Fig. 3(a). Then we consider the reso- 
nant case and tune the driving parameters to the level-crossing 
point (s()/cl>, U) = (3.832, 50). The corresponding time evo- 
lution is shown by curve 4 in Fig. 3(b), where Po(0 is always 
equal to the initial value 1, indicating the occurrence of CDT. 
But at the level-uncrossing points for non-interacting bosons 
with U - and weakly interacting bosons with U - 2, tun- 
neling still exists, as displayed by the curve 5 and curve 6 of 
Fig. 3(b). 

CDT at the avoided level-crossing points. For the avoided 
level-crossing points of partial levels, {sq/oj, U) - (2.405, 2) 
and (3.832, 52), from Eq. (14) the time evolutions of Po(0 
are plotted as the coincided curve 1 and curve 4 of Fig. 3 re- 
spectively. The invariant probability Po(t) - 1 means CDT 
of the second kind, where Eq. (10) is the superposition state 
of |i^o) and \ip2} with the same energy Eq - E2, due to 
ci - ji™ (^« ~ ") V^-^n + (^« + u)^/(^Jnkn) - for the given 

parameters. The analytical results are confirmed numerically 
from Eq. (3) for the same initial conditions and system pa- 
rameters as those of the analytical calculations. 

CDT for arbitrary values of the system parameters. If the 
system is prepared initially in a single Floquet state of Eq. 
(9), probability Pj{t) of the system in any j-th Fock state 
\i, 2 - i) is a constant for arbitrary values of the system param- 
eters in high-frequency regime, including those at the points 
with level-crossing and without level-crossing. Such an in- 
variant population means CDT of the third kind. Based on 
such CDT, we can prepared different quasi-stationary states. 
As an example, we here are interested in the initial NOON 
state |(/ro(0)) = -^[|0,2) + (-1)"+'|2,0)] of Eqs. (9) and (6). 
Under the corresponding initial conditions b'^(Q) - -b'^iQi) - 
1/ V2, Z7j(0) = 0, starting form Eqs. (11-13), time evolutions 
of Pj{t) are shown in Fig. 3(c) for cj = 50, sq/w - 2, and 
U - 2, where the values of Pq = l^oP and P2 - |CoP are 
always 1/2 and the value of Pi - |Z?()P is always 0. Thus the 
non-Floquet state of Eq. (10) actually equates the first Floquet 
state of Eq. (9), |(A(0) = |iAo(f)>, due to co = 1, ci = C2 = 0. 
Such a time-dependent state possesses the same constant pop- 
ulation with a NOON state and is called the quasi-NOON state 
thereby. For the same parameters and the initial conditions 
flo(O) = 1/ V2, fli(O) = 0, ^2(0) = -1/ numerical results 
based on Eq. (3) are shown by the solid lines of Fig. 3(c), 
which coincide with the obtained analytical results. The re- 
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FIG. 3; (Color online) Time evolutions of Fq for the parameters y - 0.5, w = 50, and 
(a) a = 2.405, V = (1,2 (coincided curve 1), f/ = 50 (curve 2), (7 = 52 (curve 3); (b) 
^ = 3.832, U = 50, 52 (coincided curve 4), (7 = (curve 5), f/ = 2 (curve 6); (c) 
^ = 2, f/ = 2 for the different initial conditions b'^(0) = -^^(O) = 1/ V2, ^',(0) = 0. 



suits consistently demonstrate that the entangled quasi-NOON 
state can be prepared by setting the initial NOON state and ap- 
plying the high-frequency driving. 

It is worth noting that Eqs. (11-13) can fit arbitrary ini- 
tial conditions of the system. Therefore, based on the differ- 
ent CDTs we can start from any initial state to construct the 
corresponding stationary-like state. As the considered exam- 
ple, for the initial conditions ^^(0) = l,b\{0) = b'^{Q) = 0, 
Eqs. (11-13) are reduced to Eqs. (14-16) which are asso- 
ciated with the initial paired state |i/'(0)) - |0, 2) and the fi- 



nal stationary-like state |i^(0) = e 



— ^mi(jjt)—i{n(ji}+u)t 



|0, 2) of Eq. 



(10). If the initial state is prepared as |i/'(0)) = |2,0) or 
|(A(0)) = (|0,2)- |2,0))/ V2, of course, from Eqs. (11-13) the 



6 



final stationary-like state reads |^fr(f)) = sm(o>tyi(ruo+u)t^2, 0> 
or the quasi-NOON state |i/r(f)) = [gl^ sin(o;o-;(«<^+«)'|o, 2) - 
sin(^o-;(«^^+«)'|2, o>]/ ^J2. These results will be used to de- 
sign the quantum tunneling switch in next subsection. 



C. Quantum tunneling switch 

From the above two subsections we know that driving field 
aff'ect tunneling dynamics dramatically, so we can design the 
quantum tunneling switch ll36ll through the high-frequency 
driving field that controls tunneling to open or close, as shown 
in Fig. 4. Let the two weakly interacting bosons occupy the 
paired state |0, 2) initially and fix the parameters co - 50 and 
U - 2. We give the scheme of quantum tunneling switch 
between different stationary-like states by adjusting the driv- 
ing strength, namely using CDT or CCT to close or open the 
tunneling. 

Switch of paired bosons tunneling. In Fig. 4(a), one can see 
that in time interval < f < fi, we let driving strength keep 
the value at the avoided level-crossing point sq - 2.405w, 
so the initially occupied state is kept, due to the CDT of the 
second kind in Fig. 3. At an arbitrarily given time t - t\, 
we change the driving strength to eq - lo) and keep this 
value until f = f3 at which the system transits completely 
from state |0, 2) to state |2,0), because of the CCT in this 
time interval. Then we return the driving strength to the initial 
value such that the CDT makes the final state the stationary- 
like state e-^sin(6;o-'('™+")'|2^0>. Thus we transport the two 
bosons from well 2 to well 1, through the quantum tunnel- 
ing switch. The switching time is just the tunneling time 
f3 - fi ~ 125(w~^) = 125/200s= 0.625s which is shown in 
the inset of Fig. 2(b). Note that in the considered case. Fig. 
2(b) points out P\(t) ~ 0. 

Switch from the initial state |2, 0) to the final quasi-NOON 
state. Before t < t2, we perform the similar operations with 
that of Fig. 4(a). At f = f2 we return the diiving strength 
to the initial value for CDT of the second kind that results in 
Po(t) ~ Piit) ~ 1 /2 for f > f2, as shown in Fig. 4(b). Here 
we have used fa as new initial time and values of Pjiti) for 
j = 0, 1,2 as the initial conditions of Eqs. (11-13) to plot 
Pj{t) for t > t2. The normalization implies P\{t2) ~ 0, be- 
cause of Poit2) + ^"2(^2) ~ 1- Therefore we have achieved the 
transition from the initial state |2, 0) to the final quasi-NOON 

state \i//it)) = eT7»in("')-<(n^^+«)'[|o,2) - e-2^™("'>|2,0)]/ V2. 
The corresponding switching time is the half tunneling time 
t2-ti s;62.5(w„i) = 0.3125s. 

In addition, making use of CDT, we can arrive at dif- 
ferent final states with the constant probabilities Po{t) - 
Pq(t), P2{t) - P2{t) for f > r, by returning the driving 
strength to the avoided level-crossing point at a different fixed 
time re (fi, t-i). Clearly, if we change the interaction strength 
from U-2toU-<x>-5Q and adjust the driving strength ei- 
ther So = 3.832<w or sq - 2co at appropiiate times, the switch- 
ing time in Fig. 4(a) can be shortened to about 6(Wq^) — 0.03s 
by CCT from the resonant interaction shown in Fig. 2(a). 




40 80 120 160 20 40 60 80 100 

t t 



FIG. 4: {Color online) Plots showing the schemes of quantum tunneling switch, 
through the time evolutions of probabilities Po(0 and ^2(0 from Eqs. (11-16) for 
oj = 50, y = 0.5. (a) In the time intervals < / < /i and t > tj, the parameter val- 
ues are taken at the avoided level-crossing point (eo/t^, U) = (2.405, 2), and in the 
time interval ti < t < tj, only the driving strength is changed to eo = 2w. (b) In the 
time intervals < / < ri and / > /2, the parameter values are taken at the same avoided 
level-crossing point with that of Fig. 4(a), and the driving strength is kept as eo - 2w in 
ti < t < li. Here r\ = 20, ti = 82.5 and = 145 are associated with the approximate 
tunnehng time - ti = 125((j^')s and its half - 'l = 125/2 = 62.5(Wo'). 



V. CONCLUSION 

We have considered two repulsive bosons confined in a pe- 
riodically driven double well and investigated the inteiplay 
between interparticle interaction and external field and its ap- 
plication to coherent control of quantum tunneling dynamics. 
Under high-frequency limit, we obtained Floquet quasiener- 
gies, Floquet states of invariant population. A full solution 
of the Schrodinger equation is constructed as a general non- 
Floquet state of slowly variant population, which is coher- 
ent superposition of the three Floquet states. The analyti- 
cal quasienergy spectra exhibit the level-crossing for reso- 
nant interaction, the avoided crossings of partial levels for 
non-resonant interaction, and the translations of level-crossing 
points caused by the photon resonance. Exploiting the co- 
herent non-Floquet states, we obtain time evolutions of the 
bosonic population, where beyond the level-crossing points, 
resonant interaction results in CCT which enhances tunneling 
rate of the paired bosons. For weak interaction, two repelling 
bosons tunnel between wells as pair and fall on a NOON state 
periodically. Three different kinds of CDT are found for the 
three-level crossing, avoided crossing of partial levels and ar- 
bitrary values of the system parameters, respectively. At the 
level-crossing points, tunneling of two bosons between wells 
is suppressed completely that means occurrence of the first 
kind CDT. Such a CDT needs larger diiving strength for a 
stronger interaction. The avoided crossing of paitial levels 
means the two-level crossing, so the superposition state of 
two Floquet states with the same energy describes CDT of 
the second kind. If the system is prepared in a single Floquet 
state initially, CDT of the third kind occurs for arbitrary values 
of the system parameters in high-frequency regime, includ- 
ing those at the level-crossing and -uncrossing points. The 
CDTs of different kinds lead to different stationary-like states 
in which the probability amplitudes are time-dependent and 
the corresponding probabilities are time-independent. 

In order to confirm the analytical results, we make direct 
numeiical simulations and demonstrate perfect agreement be- 
tween both results. As an application of these results, we 
present an interesting design scheme of quantum tunneling 
switch, through the paired-particle tunnelings between differ- 
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ent stationary-like states which contain the entangled quasi- 
NOON state, by using CDT or CCT to close or open the quan- 
tum tunneling. Such a scheme could be useful for controlling 
quantum tunneUng of two bosons in a double well. 
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